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Abstract. A lattice is called well-rounded if its minimal vectors span the 
corresponding Euclidean space. In this paper we study the similarity classes 
of well-rounded sublattices of 1? . We relate the set of all such similarity 
classes to a subset of primitive Pythagorean triples, and prove that it has the 
structure of a noncommutative infinitely generated monoid. We discuss the 
structure of a given similarity class, and define a zeta function corresponding 
to each similarity class. We relate it to Dedekind zeta of Z[i] , and investigate 
the growth of some related Dirichlet series, which reflect on the distribution 
of well-rounded lattices. We also construct a sequence of similarity classes of 
well-rounded sublattices of T? , which gives good circle packing density and 
converges to the hexagonal lattice as fast as possible with respect to a natural 
metric we define. Finally, we discuss distribution of similarity classes of well- 
rounded sublattices of Z 2 in the set of similarity classes of all well-rounded 
lattices in M 2 . 
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1. Introduction and statement of results 

Let N > 2 be an integer, and let A C M. N be a lattice of full rank. Define the 
minimum of A to be 

|A|= min ||x||, 

xGA\{0} 

where || || stands for the usual Euclidean norm on M. N . Let 

S(A) = {x e A: \\x\\ = |A|} 

be the set of minimal vectors of A. We say that A is a well-rounded lattice (abbre- 
viated WR) if 5(A) spans R^. WR lattices come up in a wide variety of different 
contexts, including discrete optimization (e.g. sphere packing, covering, and kissing 
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number problems), coding theory, and the linear Diophantine problem of Frobcnius, 
just to name a few. In particular, the classical discrete optimization problems on 
lattices can usually be reduced to WR lattices in every dimension. Distribution 
of unimodular WR lattices in M. N has been studied by C. McMullen in [17J . Also, 
the distribution of full-rank WR sublattices of Z 2 has been recently studied in [10] . 
The goal of this paper is to continue this investigation from a somewhat different 
perspective. In particular, in |10j the zeta function Cwr(z 2 )(s) of WR sublattices of 
Z 2 has been introduced, and we proved that it is analytic in the half-plane 3?(s) > 1 
with a pole of order at least two at s = 1. In Theorem 1 1.51 we establish that in fact 
the order of the pole is exactly two, where the notion of the order of the pole we use 
here is defined by (|T3|) below. To obtain this result we study the structure of the 
set of similarity classes of WR sublattices of Z 2 (Theorems 11.11 11.21 and II. 3p and 
use it to provide a simple analytic description for the Dirichlet series corresponding 
to each such similarity class (Theorem 11.41) . We then decompose Cwr(z 2 )( s ) over 
similarity classes to prove Theorem 11.51 We also discuss sphere packing density of 
similarity classes of WR sublattices of Z 2 (Theorem 11.61 and Corollary 1 1.7|) . as well 
as their distribution among all WR similarity classes in K 2 (Theorem II. 8p . 

We start out with a few words of motivation for the problems we study here. 
The similarity classes of the integer lattice Z 2 and the hexagonal lattice (defined 
in (|20p below) are very special in dimension two: these are the only two strongly 
eutactic similarity classes in K 2 , and (Ah) is the only strongly perfect similarity 
class (we define the notions of strong eutaxy and perfection at the end of Section 5, 
see (|62|) in particular; also see [16], especially chapter 16, for a detailed discussion 
of strongly eutactic and strongly perfect lattices and their properties). The distri- 
bution of sublattices of Ah is studied in [6] , while the distribution of all sublattices 
of Z 2 is well understood (see for instance [ID] , especially (52) and the beginning of 
Section 8, for a discussion of this). Studying the distribution of WR sublattices of 
these lattices is arguably even more important, since the WR property is vital in 
lattice theory. The goal of [TDJ and the current paper is to carry out this investiga- 
tion for Z 2 . We now introduce necessary notation and describe our results in more 
details. 

Recall that two lattices Ai, A2 C M. N of rank TV are said to be similar if there 
exists a matrix A in Ojv(R), the group of TV x N real orthogonal matrices, and a 
real constant a such that Ai = aAA 2 . This is an equivalence relation, which we 
will denote by writing Ai ~ A2, and the equivalence classes of lattices under this 
relation in M. N are called similarity classes. The distribution of sublattices of 1i N 
among similarity classes has been investigated by W. M. Schmidt in [21j . 

The first trivial observation we can make is that WR property is preserved under 
similarity. In other words, if two full-rank lattices Ai,A2 C are similar, say 
Ai = aAA 2 for some a €E M and A €E Ojv(M), then 

det(A 1 )-|a| A 'det(A 2 ), |A X | = |a||A 2 |, 

and Ai is WR if and only if A2 is WR. Therefore we can talk about similarity 
classes of well-rounded lattices in Mr. From now on we will write WR(il) for the 
set of all full-rank WR sublattices of a lattice tt; we will concentrate on WR(Z JV ), 
so let us write T> n and 9Jt/v for the sets of determinant and squared minima values, 
respectively, of lattices from WR(Z JV ). We will also write Cm for the set of all 
similarity classes of lattices in WR(Z W ): this is a slight abuse of notation, since 
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elements of Cn are really nonempty intersections of similarity classes of lattices in 
R N with WR(Z W ), as indicated in below when N = 2. 

In this paper we study the case N = 2. It is known that for every A £ WR(Z 2 ) 
the set S(A) has cardinality 4, and contains a minimal basis for A, which is unique 
up to ± signs and reordering (see Lemma 3.2 of [TD]). For each q £ Z>o, define 

(1) ^ = {j eQn (^' 1 ) : gcd(p,g) = i,v / 9 rr ^ez|, 

and let 

(2) S= I |J S q \ U{1}, 

where 1 is also thought of as p/q with p = q = 1. It is easy to see that the union 
in ((2J is disjoint, and each S q is a subset of the set of Farey fractions of order q in 

the interval lj ■ In Section 2 we show that the similarity classes of lattices in 

WR(Z 2 ) are in bijective correspondence with fractions p/q £ S. From now on, for 
each p/q £ S, we will write C(p, q) for the corresponding similarity class in C2, the 
set of all similarity classes of lattices in WR(Z 2 ); a formal definition of C(p,q) is 
given by ([25]) . The class C(l, 1) plays a special role: it is precisely the similarity 

class of all orthogonal well-rounded lattices, i.e. lattices of the form ^ Z 2 

for some a, b £ Z. The set C2 has interesting algebraic and combinatorial structure. 
It is not difficult to notice that the set S, which parametrizes C2, is in bijective 
correspondence with the set of primitive Pythagorean triples whose shortest leg is 
less than half of the hypothenuse. In Section 2 we explore this connection in details 
and use it to prove the following result. 

Theorem 1.1. The set C2 of similarity classes of lattices in WR(Z 2 ) has the al- 
gebraic structure of an infinitely generated free non-commutative monoid with the 
class C(l, 1) of orthogonal well-rounded lattices serving as identity. As a combina- 
torial object, C2 has the structure of a regular rooted infinite tree, where each vertex 
has infinite degree, which is precisely the Cayley digraph of this monoid. 

Remark 1.1. If G is a monoid with a generating set X, then we define its Cayley 
digraph to be a directed graph with vertices corresponding to the elements of G, 
and with a directed edge between vertices g and h if h = gx for some x £ X (see 
for instance [TH] for details and related terminology). 

We explicitly construct the monoid and the corresponding tree structure for C 2 in 
Section 2. Notice that due to Theorem 11.11 it makes sense to think of C2 as the 
moduli space of lattices in WR(Z 2 ). 



In Section 3, we discuss a more explicit parametrization of C2, which allows to see 
the structure of each similarity class C(p, q). It turns out that, although most well- 
rounded lattices are not orthogonal, all similarity classes in C2 can be parametrized 
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by a subset of lattices from (7(1, 1). More precisely, let us define a subset of Z 2 
A = | (a, 6) G Z 2 : < b < a, gcd(a, 6) = 1, 2 \ (a + 6), 

(3) and either b < a < VSb, or (2 + \/3)6 < aX, 
and consider the corresponding subset of (7(1, 1) 

(4) C"(l,l) = {(j ^ 2 eC(l,l):( a ,6)ei 

In Section 3 we prove the following theorem. 
Theorem 1.2. For eachp/q G S, there exists a unique lattice 

n=(j ffl 6 )^cui), 

where a, 6 G Z>o are given by 

(5) p = max{a 2 — 6 2 , 2a6}, and g = a 2 + b 2 , 
such that A G C(j>, q) if and only if 



(6) A = span z I x, 



yV-p 2 _p 

9 



p V q 2 -p 2 

9 9 



for some x G O. Moreover, every lattice in the set C"(l,l) parametrizes some 
similarity class C {p, q) with p, q as in in this way. 



An easy consequence of Theorem 11.21 is the existence of a lattice in each similarity 
class C(p, q) which, in a sense to be described below, generates C(p, q). First let us 
recall that given a full-rank lattice A in R 2 , its Epstein zeta function is defined by 

e a (s) = ^'iizir 2 *, 

a;GA 

where s G C, and ' indicates that the sum is taken over all x G (A/{±1}) \ {0}. 
For each such A, this Dirichlct scries is known to converge for all s with 5R(s) > 1. 
Moreover, E\(s) has analytic continuation to C except for a simple pole at s = 1. 
For more information on E\(s) and its properties see [20] . In Section 3 we also 
prove the following theorem. 

Theorem 1.3. Let G'(p,q) G £2- There exists a lattice A P)9 € C(p,q), satisfying 
the following properties. 

(1) \K p . q \ = min{|A| : A G C(p,a)} = V?, 

(%) det(A M ) = min{det(A) : A G C(p, g)} = p 7 

f5j TTie norm form of A Pi9 respect to its minimal basis is 



} P ,q( x > y) = i x% + 2a; z/vV -p 2 + ?y 2 i 



For eac/i A G C(p, g) f/iere exists U G O2QR) s^ca £/ia£ A = ^/ dct ^( A ) Jj\ p q - the 
quadratic form ( det ( A ) ) Q P: q(x, y) is therefore the norm form for A with respect to 



p 

its minimal basis, 
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(5) The Epstein zeta function of any lattice A £ C(p, q) is of the form 

and so K p . q maximizes E\(s) on C(p,q) for each real value of s > 1. 

We call Ap„ a minimal lattice of its similarity class C(p,q); it is unique up to a 
rational rotation. 

Lattices A Pj9 also determine zeta-functions of corresponding similarity classes 
C(p, q). Namely, with each C(p, q) £ C2 we can now associate two Dirichlet series, 
which incorporate information about the determinants and the minima of lattices 
in this similarity class, respectively. Specifically, define 

zt, q (s)= J2 ( det ( A )r s > z tM= E i A r 2s ' 

AeC(p,q) AeC(p,g) 

where s £ C. Our next goal is to investigate the properties of Z£ (s) and Z™ q (s) 
for each C{p 1 q) £ C2, which we do by relating them to the Epstein zeta function of 
the lattice f2 parametrizing C(p,q), as in Theorem ll.21 

We will also write (k(s) for the Dedekind zeta function of a number field K. It 
is known to be analytic for all s £ C with 3?(s) > 1 — 1/d, where d = [K : Q], except 
for a simple pole at s = 1. For more information on properties of Cr"( s ) see |15) . 
There is a standard relation between Dedekind zeta of imaginary quadratic fields 
and lattices of rank two, a special case of which we exploit here; see |23j for more 
details. 



Theorem 1.4. For each C(p,q) £ C2, 

(det(A p , 9 )) 



and 

1 1 

( 8 ) z " 9 ( s ) = ^CawW = [ApJ^ CqwW- 

We prove Theorem ll.4l in Section 3, as well. Notice in particular that C(l, 1), the 
similarity class of all lattices coming from ideals in Z[i], has Zf^s) = Z™ x (s) = 
Cnj[i](s)i since Ai,i = 1? . This fact is also discussed in [TO] , 

In [lOj we studied basic properties of the zeta- function of all well-rounded lattices 

Cwr(z 2 )(s)= E ( det ( A )r s - 

AGWR(Z 2 ) 



It also makes sense to define 



CwW*)= E i A r 2s - 

AeWR(Z 2 ) 



These two Dirichlet series carry information about the distribution of lattices in 
WR(Z 2 ) with respect to their determinant and minima values. For each similarity 
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class C(p,q) 6 C2, let us call p its determinant weight and q its minima weight. 
Theorem II .41 immediately implies that 

CwW) -EE ( det ( A ))~ s 

c( P ,q)ec 2 AeC( P , g ) 



(9) = E Oo-oowW E 3 

c( P ,q)ec 2 p ■ p/ges 



for some t/£EZ>o 



and similarly 

(10) Cwr ( z 2 )(s)=Cq W (s) E 3 



for some 



where a p is the number of similarity classes in C2 with determinant weight p, and b q 
is the number of similarity classes in C2 with minima weight q; notice that b q = |«S 9 |, 
where S q is as in JT]). In fact, let us write 

(id w<(-)= E i= E ^ 

C( P ,g)eC 2 y P ■. p/qes 

for some g£Z>o 

and 

(12) w TO (-)= E i= E 

C(p,g)GC 2 y 9 : p/geS 

for some p£Z>o 

We will call Wd(s) and W m {s) determinant and minima weight enumerators, re- 
spectively. Therefore the question of distribution of lattices in WR(Z 2 ) is linked to 
understanding the basic analytic properties of Wd{s) and W m {s). In Section 4 we 
use an approach different from that of [lOj to prove the following result. 

Theorem 1.5. Let the notation be as above, then Wd{s) and W m (s) both have 
simple poles at s = 1 and are analytic for all s € C with 5R(s) > 1. Therefore 
Cwr(z 2 )( s ) an d Cwr(z 2 )( s ) have poles of order two at s = 1 and are analytic 

for all s eC with 3?(s) > 1. 

We should point out that we are using the notion of a pole here not in a sense that 
would imply the existence of an analytic continuation, but only to reflect on the 
growth of the coefficients. More precisely, for a Dirichlct series Y^Li c nn~~ s , we say 
that it has a pole of order ji at s = So, where fi and sq are positive real numbers, if 

00 

(13) 0< lim \s- s \»^>2\c n n- s \ < 00. 

s ^ s o n=l 

Notice that Theorem 11.51 in particular improves slightly on the result of Theorem 
1.5 of [10]. The approach we use in Section 4 to prove Theorem 11.51 uses bounds 
on coefficients of weight enumerators Wd{s) and W m (s) by coefficients of Dirichlct 
series associated with the set of primitive Pythagorean triples, which have Euler 
product expansions. 



A standard object of lattice theory is a sphere packing associated with a lattice, 
and a classical problem is to determine the optimal packing density among lattices 
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in a given dimension (see [9]). This problem has been solved in dimension two; in 
fact, it is not difficult to show that maximization of packing density can be restricted 
to WR lattices. Here we will discuss the circle packing density corresponding to 
lattices in WR(Z 2 ), investigating how "close" can one come to the optimal packing 
density in dimension two with such lattices. For these purposes, let us write (A) 
for the similarity class of any lattice A in R 2 , so that 

(14) (A) = {aUA : a e R> , U G 2 (R)}- 
Then for each p/q E S, 

(15) C(p, 9 ) = (A M )nWR(Z 2 ). 
For a lattice A in R 2 define 

C / \ X t y\ \ 

(16) #(A) = min < arcsin - — — — - ) : x, y is a shortest basis for A 

v ' K ' I \\\x\\\\y\\) 

By a shortest basis x,y of A we mean here that a; is a minimal vector of A, and 
y is a vector of smallest Euclidean norm such that x, y is a basis for A. By a well 
known lemma of Gauss, 9(A) S [§,§] (see [M]). It is easy to notice that 9(A) 
remains constant on (A), so we can also write 9((A)). If x,y is a shortest basis for 
A with the angle between x and y equal to 9(A), then 

det(A) = ||as||||2/|| sin (9 (A), 
and so if A is well-rounded, then = \\y\\ = |A|, and so 

(17) det(A) = |A| 2 sin 0(A). 

It is easy to see that two well-rounded lattices Aj , A 2 C R 2 are similar if and only 
if 0(Ai) = 9(A 2 ), i.e. if and only if 



sin^(Ai) = sin0(A 2 ) € 



2 ' 



and so similarity classes of well-rounded lattices in R 2 are indexed by real numbers 



in the interval 



2 -1 



Let Sim(R 2 ) be the set of all similarity classes of well- 
rounded lattices in R 2 , and for every two (Ai) , (A 2 ) G Sim(R 2 ) define 
(18) d.(Ai,A a ) - \sm9(A 1 )-sm9(A 2 )\. 

It is easy to see that d s is a metric on Sim(R 2 ). If A is a well-rounded lattice in 
R 2 , then the density of circle packing given by A is 

tt|A| 2 tt 



(19) 5(A) = 4det(A)-4sin0(A)' 

by (|17p , and so it depends not on the particular lattice A, but on its similarity class 
(A). Moreover, (TIT))) implies that the smaller is sin 9(A) the bigger is S(A). Indeed, 
it is a well known fact that the similarity class (A^) gives the optimal circle packing 
in dimension two, where 



l 

2 \ ^2 



(20) A h =Q ^jZ 

is the two-dimensional hexagonal lattice, and sin9(Ah) = The lattice A^ also 
has the largest minimum among all lattices in R 2 with the same determinant, and 
minimizes Epstein zeta function for all real values of s > 1 (see [7])- However, 
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(Ah) H WR(Z 2 ) = 0. How well, with respect to the metric d s on Sim(R 2 ), can we 
approximate the similarity class (Ah) with similarity classes of the form (A p , q ), i.e. 
with similarity classes that have a nonempty intersection with the set WR(Z 2 )? 
This question is especially interesting since, in contrast to the two-dimensional 
situation, the three-dimensional counterpart of Ah, the face-centered cubic (fee) 
lattice which maximizes sphere packing density in R 3 , is in WR(Z 3 ). Our next 
result addresses this question. 

Theorem 1.6. There exists an infinite sequence of similarity classes (A Pk , qk ) such 
that 

( A P fc ,?J — ► ( A h) , as k -> oo, 

with respect to the metric d s on Sim(R 2 ). The rate of this convergence can be 
expressed by 

(21) - I < d s (A h ,A Ph , qk ) < ' 



where q k = 0(14 ) as k — > oo. Moreover, the inequality V21)) is sharp in the sense 
that 

(22) -^=-<d s (A h ,A p , q ), 

for every similarity class of the form (A p>q ) ^ (Ai.i). For the similarity class of 
orthogonal well-rounded lattices (A^i) = (j^ 2 ) , we clearly have d s (A/,, Z 2 ) = • 

Corollary 1.7. Each similarity class (A Pk ^ qk ) of Theorem ] 1 . 6\ gives circle packing 
density 5 Pk . qk such that 

(23) 6 (An) ( 1 | ) <S Ph , gk <S(Ah)( - \ 92 ) , 

\ ' 723x(13.928)''- 1 / V 723x(13.947) fc - 1 / 

where 5(Ah) = -^7= = 0.9069... is the circle packing density of Ah- 

We prove Theorem 11.61 and Corollary 1 1.71 in Section 5. Notice that a well-rounded 
lattice in M 2 has a rational basis, i.e. a basis consisting of vectors with rational 
coordinates, if and only if it belongs to a similarity class (A pq ) for some p, q. 
Therefore results of Theorem ll.6l and Corollary QT7] can be interpreted as statements 
on best approximation to A^ (and hence best circle packing) by well-rounded lattices 
in M 2 with rational bases. As we will see in Section 5, this just comes down to finding 
best approximations to by fractions | where (p, \J q 2 — p 2 , q) is a primitive 

Pythagorean triple with yj q 2 — p 2 < q/2. In fact, a similar approximation result 
holds for all WR lattices in R 2 , not just Ah- 

Theorem 1.8. The similarity classes of WR sublattices of I 2 are dense in the set of 
all similarity classes of WR lattices in R 2 , in other words the set {(A p-9 ) : p/q € S} 
is dense in Sim(R 2 ) with respect to the metric d s . Moreover, for every A 6 Sim(P 2 1 

2V2 



there exist infinitely many non-similar lattices A p q S WR(Z 2 ) such that 



(24) d.(A,A M )< 
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We derive Theorem 11.81 in Section 6 as an easy corollary of a theorem of Hlawka 
on Diophantine approximation with quotients of Pythagorean triples, and discuss 
equidistribution of {(A p . q ) : p/q G S} in Sim(R 2 ). As a side remark in Section 6, 
we also use Hlawka's result to approximate points on a rational ellipse by rational 
points on the same ellipse. Notice that Theorem 11.81 does not include Theorem 11.61 
as a special case, since the approximating constants in Theorem [L6] are sharper and 
the proof is constructive unlike that of Theorem 1 1.81 We are now ready to proceed. 



2. Parametrization by Pythagorean triples 

Notice that if a lattice A G WR(Z 2 ), then cos 9(A), sin 9(A) G Q >0 , where 9(A) is 
defined in (fTB)) . and therefore we can index similarity classes of lattices in WR(Z 2 ) 
by fractions p/q G S, where S is as in ([2]), so for each such p/q the corresponding 
similarity class C(p, q) G Ci is a set of the form 

(25) C(p,q) = | A 6 WR(Z 2 ) : sin 0(A) = J 



For each p/q G S, define t — y/ 'q 2 — p 2 G Z. Then it is easy to notice that 

q \fiq 
0<t<l<^<p<q, 

and t 2 + p 2 = q 2 with gcd(t,p, q) — 1. In other words, the set 5, and therefore 
the set C2 of similarity classes of lattices in WR(Z 2 ), is in bijective correspondence 
with the set of primitive Pythagorean triples with the shortest leg being less than 
half of the hypothenuse. 
Let 

ty = {(x,y,z) ■ x,y,zeZ >0 , 2\y, gcd(x, y, z) = 1, x 2 + y 2 = z 2 } 

be the set of all primitive Pythagorean triples, and let 

V = {(x,y,z)e^ : min{z,y} < z/2} U{ (1,0,1)}. 

Notice that we include (1,0, 1) in V '. although it is traditionally not included in ^. 
Then p/q € S if and only if either (t,p,q) € V or (p,t,q) G V . In other words, 
elements of V can be used to enumerate similarity classes of lattices in WR(Z 2 ). We 
will use this approach to provide a convenient combinatorial description of elements 
of C2. Define matrices 

/l -2 2\ /l 2 2\ /-l 2 2\ 

(26) A=[2 -1 2 ,5= 2 1 2 , C = -2 1 2 e GL 3 (Z), 

\2 -2 3/ \2 2 3/ \-2 2 3/ 

and let G = A, B,C) be the non-commutative monoid generated by A,B,C 
with the 3x3 identity matrix I3. Let us think of elements of *p as vectors in Z 3 , 
and for each MgG define the corresponding linear transformations 

(27) M(x,y,z) = M ly 

It is a well known fact that for every (x, y, z) G *p, A(x, y, z), B(x, y, z), C(x, y, z) G 
^p. Moreover, every (x, y, z) G ^3 can be obtained in a unique way by applying 
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a sequence of linear transformations A,B,C to (3,4,5), the smallest triple in 
(this construction is attributed to Barning [4]; also see pQ, [18] ). This means that 
(|27p defines a free action of G on the set of primitive Pythagorean triples by left 
multiplication. The set has the structure of an infinite rooted ternary tree with 
respect to this action, as described in pQ ; this in particular implies that G is a free 
monoid. In fact, this tree (see Figure 1 below) is precisely the Cayley digraph of G 
with respect to the generating set {A, B, C}. 



(3,4,5) 




(5,12,13) (21,20,29) (15,8,17) 




Figure 1. Ternary tree representation for *}3 

We can extend this construction by considering the set = *J5 U {(1,0,1)} 
(compare with [2\). It is easy to notice that A(1,Q, 1) = 5(1,0,1) = (3,4,5) and 
C(1,0, 1) = (1,0, 1), and hence every (x,y,z) 6 can be obtained by applying a 
sequence of linear transformations A, B, C to (1, 0, 1). Such a sequence is no longer 
unique, hence action of G does not extend to however a shortest such sequence 
is unique up to multiplication on the left by either BA~ X or AB^ 1 . 

Notice that V c Let 

H = Ia 2 NB, A k B, ABNB, C 2 NB, C k B, CBNB 

{AC) k A 2 NB, (ACfABNB, {AC) k AB, 

(28) {CA) k C 2 NB,{CA) k CBNB,{CA) k CB : N e G, k e Z >0 }. 

It is clear that H is a subsemigroup and H' = H U {J3} is a submonoid of G. Let 
us define the image of under G to be 

GW = {M(x,y,z) : M e G, (x, y, z) g 

and similarly for the images CP, -ff^', and HV. 

Lemma 2.1. = HV = P \ {(1, 0, 1)}. 

Proof. First we will prove that H^' C P. It is clear that (1,0,1) i Hty. Let 
M € if, then there exists some N £ G such that one of the following is true: 

(1) M = A 2 N, 

(2) M = ABN, 

(3) M = C 2 7V, 

(4) M = CBN , 

(5) M = {AC) k A 2 N, where fc G Z >0 , 
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(6) M = (AC) k ABN, where k e Z >0 , 

(7) M = (CA) k C 2 N, where k G Z >0 , 

(8) M = (CA) fe CBiV, where k e Z >0 . 

Let (a;, y, z) e and write (a;', y', z') = iV(a;, y, z), where N is as above. Then, in 
case (1) 

I x' - Ay' + 4z' \ 
M(a:, y, z) = A 2 (x', y', z') = 4a;' - 7y' + 8z' , 

\4a;' - 8y' + 9z' ) 

where 

i(4x' - 8y' + 9z') = 2x' - Ay 1 + ^z' > x' - Ay 1 + Az' , 
hence M(x,y, z) e V . In case (2) 

/ x' + Ay' + Az' \ 
M(x, y, z) = AB{x', y' , z') = 4a;' + 7y' + 8z> , 

\4a;' + 8?/ + 9z'J 

where 

1 9 

-(4a;' + 8y' + 9z') = 2x' + Ay 1 + -z' > x' + Ay 1 + Az', 

hence M(x,y, z) e V . In case (3) 

/-7a;' + 4y' + 8z'\ 
M(x 1 y,z) = C 2 (x' 1 y' 1 z') = -Axf + y' + Az' , 

\-8x' + Ay' + 9z'J 

where 

i(-8x' + Ay' + 9z') = -Ax' + 2y' + ^z' > -Ax' + y' + Az' , 
hence M(x,y, z) e V . In case (4) 

/7a;' + 4y' + 8z'\ 
M(x,y,z) = CB(x',y',z')= Ax' + y' + Az' , 

\8x' + Ay' + 9z'j 

where 

1 9 

-(8a;' + Ay' + 9z') = Ax' + 2y' + -z' > Ax' + y' + Az' , 

hence M(x, y, z) e V . 

For cases (5) and (6), let 

(x2,y2,z 2 ) = M(x,y,z) = (,4C) fe (a;i,yi,zi), 

where 

fx'- Ay' + Az' \ 
(x u y u z 1 )=A 2 (x',y',z') = Ax' - 7y' + 8z' 

\Ax' - 8y' + 9z'J 

in case (5), and 

/ x' + Ay' + Az' \ 
(x 1 ,y 1 ,z 1 ) = AB(x',y',z') = Ax' + 7y' + 8z' 

\Ax' + 8y' + 9z'J 

in case (6). It is not difficult to notice that Xi = min{a;2,y2}, and 

Z2 . fZl \ 

y = X2 + {y- X1 )- 
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Therefore (x2, 2/2, Z2) € V if and only if xi < zi/2, which is true in both cases, (5) 
and (6). On the other hand, 

/ x' + 2y' + 2z' \ 
B(x',y',z')= 2x' + y' + 2z' , 
\2x' + 2y' + 3z'J 

and 

/ -x' + 2y' + 2z' \ 
C(x',y',z')= -2x' + y' + 2z' , 
\-2x' + 2y' + 3z'J 

which implies that (AC) k B(x', y', z'), (AC) k C(x', y', z') £ V for any (x' , y' , z'). 
For cases (7) and (8), let 

(x 2 ,y2,z 2 ) = M(x,y,z) = {CA) k (x 1 ,y 1 ,z 1 ), 

where 

/-7x' + Ay' + 8z'\ 
(x 1 ,y 1 ,z 1 ) = C 2 (x',y',z')= -Ax' + y' + Az' 

\-8x' + Ay' + 9z'J 

in case (7), and 

hx' + Ay' + 8z'\ 
(x 1 ,y 1 ,z 1 )=CB(x',y',z , )= Ax' + y' + Az' 

\8x' + Ay' + 9z'J 

in case (8). It is not difficult to notice that j/2 = mining, 2/2}, and 

z 2 (z 1 \ 

Therefore (X2, j/2, Z2) € V if and only if y\ < Zi/2, which is true in both cases, (7) 
and (8). On the other hand, 

/ x' + 2y' + 2z' \ 
B(x',y',z')= 2x' + y' + 2z' , 
\2x' + 2y' + 3z'J 

and 

fx'- 2y' + 2z' \ 
A(x',y',z')= 2x'-y' + 2z' , 
\2x' - 2y' + 3z'J 

which implies that (CA) k B(x',y',z'), (CA) k A(x',y',z') £ V for any (x 1 , y' , z'). We 
have shown that HT C H<$' CV\ {(1, 0, 1)}. 

To finish the proof of the lemma, we will show that V \ {(1, 0, 1)} C HT. Notice 
that it is in fact sufficient to show that for each (x, y, z) 6 P\ {(1, 0, 1)} there exists 
M E H such that (x,y, z) = M(1,0,1). We know that there exists N £ G such that 
(x, y, z) — N{3, 4, 5), and so (x, y, z) — NB(1,0, 1). First notice that N cannot be 
of the form BN' for some N' e G. Indeed, suppose it is, then 

fx' + 2y' + 2z' \ 
{x,y,z) = B(x',y',z')= 2x' + y' + 2z' , 

\2x' + 2y' + 3z'J 

where (x', y', z 1 ) = N'(x, y, z) e but 

i(2x' + 2y' + 3z') = x' + y' + ^z' < min{x' + 2y' + 2z', 2x' + y' + 2z'}, 
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which contradicts the fact that (x, y, z) £ V . Similarly, from the arguments above 
it follows that N cannot be of the form (AC) k BN', {AC) k CN', (CA) k BN', or 
{CA) k AN' . The only options left are those described in cases (1) - (8) above, 
which means that M = NB £ H . Therefore V \ {(1, 0, 1)} C HV, which completes 
the proof. □ 

Theorem 2.2. H' is a free infinitely generated monoid, which acts freely on the set 

V by left multiplication. With respect to this action, V has the structure of a regular 
rooted infinite tree, where each vertex has infinite degree (see Figure 2 below); this 
is precisely the Cayley digraph of H' . 

Proof. H' is a submonoid of G, which is a free monoid, hence H' must also be free 
by the Nielsen-Schreier theorem (see for instance [19]). To see that H' is infinitely 
generated, consider for instance the set {AB k : k £ Z>o} of elements of H' . Since 
A, B k ^ H' for any k £ Z>o, it is clear that no finite subset of H' can generate all of 
the elements of the form AB k : if this was possible, there would have to be relations 
between elements of H' , contradicting the fact that it is free. Therefore H' must 
be infinitely generated, and so its Cayley digraph is a regular rooted infinite tree, 
where each vertex has infinite degree, and the root corresponds to I3. 

By Lemma 12.11 we know that H'V = V . Moreover, we know that for each 
(x,y,z) 6 V there exists a unique element N 6 G such that iV(3,4, 5) = (x,y,z), 
and hence NB is the unique element in H such that NB(1,0, 1) = (x,y,z). This 
means that H' acts freely on V. Then we can identify (1,0,1) € V with ^3 £ 
H', and each (x,y,z) £ V with the corresponding unique NB £ H' such that 
NB(1, 0, 1) = (x, y, z), which means that with respect to the action of H' the set 

V has the structure of the Cayley digraph of H' with respect to an appropriate 
generating set. □ 




Figure 2. Infinite-degree tree representation for V 



Corollary 2.3. The set C2 of similarity classes of lattices in WR(Z 2 ) has the 
structure of a non- commutative free infinitely generated monoid. Specifically, it is 
isomorphic to H' . 

Proof. We will identify C2 with H' in the following way. From Theorem 12.21 we 
know that there exists a bijection ip : V — > H' , given by 

ip(x, y, z) = M, such that M(l, 0, 1) = (x, y, z), 
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for each (x, y,z) E V with ip 1 : H' — > V defined by 

(p~ 1 (M)=M (1,0,1), 

for each M G if'. 

On the other hand, there also exists a bijection ip : C 2 —> V, given by 



^(C(p,g)) 



Wq 2 -p 2 ,p , q) if2|p 

{p^q 2 -p 2 ,q) if 2fp, 
for each C(p, q) S C2 with : P — > C2 defined by 

i/j^ 1 (x,y, z) = C(p,q), where p = max{x, y}, q = z, 
for each (x,y,z) E V. Therefore we have bijections ^ : C2 — * H' and ((/??A) _1 = 

We can now define a binary operation * on C 2 as follows: for every C(pi, q±) and 
C(p 2 ,q 2 ) in C 2 , let 

(29) C(pi, gi) * C( P2 , q 2 ) = Tjr 1 (< P Tp(C( Pl ,q 1 ))<piP(C(p 2 ,q 2 ))(l,0, 1)) . 

It is easy to see that C 2 is a free non-commutative monoid with respect to *, which 
is isomorphic to H' via the monoid isomorphism (pip : C 2 — > H', and ((pip) -1 = 
C(l, 1) e C2 is the identity. Hence the tree in Figure 2 is the Cayley digraph of C 2 
with respect to an appropriate generating set. This completes the proof. □ 

Now Theorem 11.11 follows by combining Theorem 12.21 with Corollary |2.3l 



3. Similarity classes and corresponding zeta functions 

In this section we discuss the structure of similarity classes C(p,q), as well as 
the properties of associated zeta functions. Our first goal is to prove Theorem 11.21 
For each p/q E S, define 

M 2 (p 7 q) = {xe (Z 2 /{±1}) \{0} : x iy /q2-p> = x 2 p (modq), 

x 2 \J q 2 — p 2 = —x\p (mod<7)|. 

Notice that x E Wl 2 (p, q) if and only if 

r /VZEZ _p \ ] 

(30) A(x) := span z i x, he U C(p, q). 

Hence lattices in the similarity class C(p,q) are in bijective correspondence with 
points in 9Jl 2 (p 1 q). 

Lemma 3.1. Let p/q E S. The congruence relations 

(31) X\ \/ q 2 — p 2 = x 2 p (mod q) 
and 



(32) x 2 \/ q 2 — p 2 = — x\p (mod q) 

are equivalent, meaning that 

m 2 (p,q) = {xe (Z 2 /{±1})\{0} : x^q 2 -p 2 = x 2 p (modg)}. 



ON SIMILARITY CLASSES OF WELL-ROUNDED SUBLATTICES OF Z 2 



15 



Proof. Recall that gcd(p, q) = 1. Also q 2 — p 2 = (p — q)(p + q), and 
gcd(g -p,q) = gcd(g +p,q) = gcd(p, g) = 1, 



therefore gcd(g 2 — p 2 ,q) = gcd(-\/g 2 — p 2 , g) = 1. Hence 

xi v / g 2 — p 2 = X2P (mod g) 

if and only if 

X2P\/ q 2 — P 2 = a?i (g 2 - p 2 ) = —xip 2 (modg), 



which happens if and only if —x\p = X2 \/ q 2 — p 2 (modg). □ 

For each p/q G S, define c(p, g) to be the unique integer such that < c(p, q) < 
q — 1 and 

(33) c(p, g)p = \] q 2 - p 2 (mod g). 

Then, by Lemma |3.1[ for every x G 9^2 (p, g) we have X2 = c(p, g)xi (modg), 
meaning that £2 = c(p, q)x\ + yq for some y G Z. In other words, 3Jt2(p, g) can be 
presented as 

(34) Wfe^-IU^ + w) : (;)e(*/{±l})\{0}}. 
Define 



1 0\ _ 2 



(35) = q) Z > 
so that m 2 (p,q) = (n(p,g)/{±l})\{0}. 

Lemma 3.2. f2(p,g) G (7(1,1), and |£7(p,g)| 2 = det(fi(p,g)) = g. In fact, each 

m i/je sei C"(l, 1) is 0/ i/ie /orm f2(p, g) /or p, g satisfying |5p. 

Proof. First fix a lattice Q(p, g). It is a well known fact that there exist unique 
relatively prime a > b G Z>o of different parity such that either p = a 2 — b 2 
or p — 2ab, and g = a 2 + 6 2 (this is the standard parametrization of primitive 
Pythagorean triples, see for instance [2Zj)- The fact that ^g < p < q ensures that 
(a, b) G A. Then the lattice 

°-(: ~'y 

is in C"(l, 1), and |r2| 2 = det(f2) = g. We will now show that f2(p, g) = 17. Since 
gcd(a,6) = 1, there exist 31,32 G Z such that 

(36) gia-g 2 b=l. 

Let 7 = 316 + 32a, and notice that 

'a -b\ f 9l b\ = h ^ 
b a J \g 2 a) \j q, 

and det ( ? | = det y \ = q. so 
\b a J \j q) 

(37) n=(i °> 
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Notice that + 1 is divisible by q. Indeed, ([36]) implies that g\ = 92 ^ +1 , and so 

b 2 

7 2 + l = (g 1 b + g 2 a) 2 + l = q 2 a 2 + 2g 2 b(g 2 b+l) + ^r(g 2 b+l) 2 + l 

or 

/ 529 + 2g 2 b+l\ 2 2 

(38) = 9 I ^ )=9(5i+.9 2 )- 

Moreover, we can ensure that < 7 < g — lby replacing 7 with 7 + gm for some 
m G Z, if necessary: it is easy to see that 7 2 + 1 will still be divisible by q, and (f3"T)) 
will still hold. We will now show that 7 = c(p, q). Notice that 

(39) F y , q (x, y) = i 1 ^-] x 2 + 2jxy + qy 2 



q 

is an integral binary quadratic form with discriminant -4, hence it is equivalent to 

G(x,y) = x 2 +y 2 , 
since the class number of -4 is one. In fact, it is easy to verify that 

F -rA x , V) = G(gix - by, g 2 x - ay), G(x, y) = F ltq {ax - by, g 2 x - giy). 



Let us write t = a 2 — b 2 , so either p — t or p = y 1 q 2 — t 2 , then 

q = G(a, b) = F 7>9 (t, k) = f 1 +1 ) t 2 + 2-ftk + qk 2 , 



where k = g 2 a — g\b. Therefore 

7 2 i 2 + 2 7 t/c + q 2 k 2 = q 2 - t 2 , 

meaning that 

(40) 7 2 i 2 =q 2 -t 2 (mod q). 

Notice that gcd(7, q) = 1, since 9K7 2 + 1), so gcd(7, q) must divide 1. Therefore, if 
p = t, then (|40|) implies that 

(41) jp = \/ q 2 — p 2 (modg). 

If, on the other hand, p = \J q 2 — t 2 , then (|4U|) implies that 71/ q 2 — p 2 = p (mod q), 
meaning that 7 2 y/q 2 — p 2 = jp (mod q), but on the other hand 7 2 = — 1 (modg), 
and so 

— \J q 2 — p 2 = -fp (mod q) . 

By Lemma [3.11 this last congruence is equivalent to (l4Tj) . We conclude that < 
7 < q — 1, and 7 satisfies (|4ip. which means that 7 = c(p, q), and so f2 = £l(p, q). 

In the opposite direction, assume that 

and define q — a 2 + b 2 , p = max{a 2 — b 2 ,2ab}. Then the fact that (a, b) G A 
ensures that ^q < p < q, i.e. p/q G S. It is not difficult to notice that for p/q G S, 
p = a 2 — b 2 if and only if a > (2 + \/3)b, and p — 2ab if and only if b < a < \/3b. The 
argument identical to the one above now shows that Q — £l(j>,q). This completes 
the proof. □ 
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Remark 3.1. It is not difficult to conclude from an argument very similar to the one 
in the proof of Lemma 13.21 that all binary integral quadratic forms of discriminant 
-4 are of the form Ky iq (x, y) as in Q39p for some odd positive integer q which is not 
divisible by any prime of the form 4fc + 3 and an integer 7 (positive or negative) 
such that 7 2 + 1 is divisible by q. This statement is essentially equivalent to the fact 
that there is a bijection between ideals of the form (da + dbi) in Z[i] with a > b > 0, 
d > 0, and Pythagorean triples (d 2 (a 2 - b 2 ), 2d 2 ab, d 2 (a 2 + b 2 )). 

Proof of Theorem ! 1 . 2[ Fix a similarity class C(p,q) £ C2 for some p/q £ S. For 

each A 6 C(p, q), we have A = A(x) as defined by (1301) for some x = \ , s ] 

\c(p 7 q)x + qyj 

where x,y € Z, hence £e G f2(p, g). On the other hand, for each x £ Cl(p,q), 
the corresponding lattice A(x) is easily seen to be in C(p,q). Combining this 
observation with Lemma 13.21 completes the proof of the theorem. □ 



Proof of Theorem \1.3[ We first introduce the notion of a minimal lattice in each 
similarity class C(p, q). Let x(p, q) £ Q(p, q) be such that \\x(p, q)\\ = \Q(p, q)\, and 
let A(x(p, q)) be defined by ([30]) : we will call this lattice a minimal lattice of the 
similarity class C(p, q) and will denote it by A Pig . By Lemma |3.2[ we have 

(42) |fi(p,?)| = ||a:(p,«)||=Vg. 

On the other hand, since -s/3/2 < p/q < 1, meaning that the angle between x(p, q) 
and the other minimal basis vector given in (|3"0|) is arcsin(p/<7) £ 7r/2), a 

well known lemma of Gauss (see [3] or [TU]) implies that |A Pj9 | = \\x(p, q)\\, and so 
l^p.gl = %/*?■ Therefore 

(43) det(A p , g ) = |A P:9 | 2 V - = p = min{det(A) : A £ C{p, <?)}. 

Moreover, a straight-forward computation shows that the norm form of A Pj9 with 
respect to its minimal basis is 



2 1 qy 2 - 



(44) Q p , q (x,y) = (x,y) ( y^T^ ^ ^) Q = 9^ + 2^^?^ 

Notice that the minimal lattice of a similarity class may not in general be unique, 
however it is unique up to a rational rotation, and so for our purposes it suffices to 
pick any one of them. 

Next, let A £ C(p,q), then A ~ A Pi9 , and so there must exist a £ R>o and 
U £ 2 (K) such that A = aUA p , q . Then det(A) = a 2 p, and so a = y^MAl > 1. If 
we write A and and A p q for the minimal basis matrices of A and A p q respectively, 
then A = / dct ( A ) jj A PtQ , and the norm form of A with respect to this minimal basis 



is 



det(A) t fx\ det(A) 



Q A (x,y) = {x,y)A t A [ = -^-^(x 1 y)A t pq A p , q = —±-LQ M (x lV ). 
Epstein zeta function of A is therefore given by 

E A ( S ) = J2' \\*\r 2s = £' QA(x,yr s = (^£^Y E' Q P A^y)- s 
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Then (j43|) implies that for every fixed real value of s > f, E\(s) achieves its 
maximum on C(p,q) when A = A Pj9 , and it does not achieve a minimum since 
there exist lattices in C(p,q) with arbitrarily large determinants. This completes 
the proof of the theorem, rj 



Proof of Theorem \l-4\ We now derive the properties of the Dirichlet series corre- 
sponding to each C(p, q). Fix a similarity class C(p, q) e Ci- By Theorem II. 2 [ each 
A e C(p, q) is of the form A(cc) for some x E f2(p, q). As in the proof of Theorem 
1.31 above, a well known lemma of Gauss (see [3] or [10]) implies that |A(a;)| = 
Since also, by Lemma 3.2 of [TU], the set of minimal vectors of A(x) is precisely 



V q 2 -p 2 



±X,±\ q I X 

p y<r-p 2 
q i 

it follows that A(a?i) = A{x^) if and only if X\ = ±X2- Therefore 

Z p^ s ) = E l A l^= E' NI~ 2s = £n (M) (s). 

AeC(p,q) xen(p,q) 

Now Theorem 1 1 . 2 1 readily implies that there exists U € such that 

n(p, ? ) = Er(f ^)z», 

which means that Sq( Pi9 ) (s) is equal to the Epstein zeta function of ( ) J 



Hence 



( S ) = £n (p , g) (s) = i E' INI" 3 ' = i<Q(0(«). 



which proves ([5]) . Now recall that for each A G C(p, q), 

det(A) = |A| 2 sin0(A) = |A| 2 ^. 

Q 

Then © follows. □ 



Finally, we present a simple lemma, which is not related to the rest of this section, 
but is of some independent interest. Recall that the dual of a lattice A in M. N is the 
lattice A* , defined by 

A* = {x e M. N : V y G A, x l y E Z}. 

Lemma 3.3. Let C(p,q) E C2, and suppose that A E C(p,q). Then det(A)A* S 
C(p,q). 

Proof. Let x, y be the minimal basis for A, and write A — \ Xl ^ 1 ) for the 
corresponding basis matrix. Then 

(«) (A- 1 )* = j — TTT ( V2 ~ X2 ) 

y ' y ' det(A) \-yi x x J 
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det(A)A* = (J Q ^ A, 



and so A - A*, since G 2 (M). □ 



4. Weight enumerators Wd(s) and W m {s) 



In this section we will discuss in more details some properties of the Dirichlet 
series Wd(s) and W m (s) as defined in l|lip and (|12p. respectively, and will prove 
Theorem 11.51 Recall that we write a p and b q for the coefficients of Wd(s) and 
W^s) respectively as defined in Section 1. The following formulas for a p and b q 
are immediate from Theorem 11.21 and the definition of the set A in © . 

Lemma 4.1. For each p such that p/q € 5 /or some q G 1>o, 



Notice that the expression for a p in Lemma |4. II is similar in spirit to the function 
(3 defined in [10], in particular it can also be bounded in terms of Hooley's A- 
function. On the other hand, we can obtain simple explicit bounds for a p and b q 
from our Pythagorean tree construction in Section 3. For each p, q € Z>o, define 
L(p) to be the number of primitive Pythagorean triples with a leg p, and H(q) to 
be the number of primitive Pythagorean triples with the hypotenuse q. There are 
well known formulas for L(p) and H(q) (see [5], p. 116): if p, q > 1, then 



b q < H(q) when p, q > 1, and a± = &i = 1. One can ask how good are these 
bounds? We will now show that the correct order of magnitude of the bound for 
both, a p and b q , in the sense that the corresponding Dirichlet series has the same 
behavior at s = 1 as Wd{s) and W m (s), is given by H and not by L. Namely, define 



a p = \{(m, n) G A : p = max{m 2 — n 2 , 2mn}} 
and for each q such that p/q £ S for some p G Z>o, 

b q = \{(m, n) € A: q = m 2 + n 2 }\ . 
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Lemma 4.2. H(s) has a simple pole at s = 1 and is analytic for all s € C with 
3?(s) > 1. Moreover when SR(s) > 1, H(s) /ias an Euler product type expansion 



(46) 



ww=| n 



p + 1 



Z=l(mod4) 

where the product is over primes I . 
Proof. Let us define 

V\ = {n € Z>o : n is only divisible by primes which are = l(mod4)} 
Then notice that, as in the proof of Lemma 8.1 of [TO] . 

2«(«) 



2H(s) = £ 



n E 2 ^" = n h+ 2 E^ 

Z=l(mod4) \fc=0 / ;=l(mod4) \ k=l 

2 A -r-r l S + 1 



— ks 



n 



/=l(mod4) 



1 - ?- 



1 



n 



i=l(mod4) 



whenever this product is convergent, where I is always prime. The fact that 7i(s) 
has a simple pole at s = 1 and is analytic for all s G C with 5R(s) > 1 then follows 
immediately from Lemma 8.1 of [10 . □ 

Proof of Theorem \1.5[ First of all notice that since -^q < p < q, we have 




E 

c( P ,9)ec 2 



C(p,9)ec 2 



* E 

C(p, 9 )ec 2 



meaning that Wd{s) and f m (s) must have poles of the same order and the same 
half-plane of convergence. Since b q < H(q), Lemma |4~21 implies that W m (s) has at 
most a simple pole at s = 1, and is analytic when SR(s) > 1. On the other hand, 
Theorem 1.5 of [10] implies that Cwr(z 2 )(s) has at least a pole of order two at 
s = 1, meaning that, by (J9j) , Wd(s) must have at most a simple pole at s = 1. This 
means that both, Wd(s) and lf m (s), have simple poles at s = 1 and are analytic 
when 3?(s) > 1, and therefore, by ([9|) and fTO]) . Cwr(z 2 )(s) an d Cwr(z 2 )( s ) both 
have poles of order two at s = 1 and are analytic when 5R(s) > 1. This completes 
the proof. □ 



Remark 4.1. Notice that Theorem 11.51 combined with Lemma I4T21 implies that the 

Dirichlet series E( x , y , z)e v max {x, y f and £( w )e<p; SS^W have poles of the 
same order 1 at s = 1. This fact could be roughly interpreted to mean that the 
sets V and ty' are comparable in size, i.e. that "most" primitive Pythagorean 
triples correspond to similarity classes of lattices from WR(Z 2 ). In other words, 
the imposed condition that the shortest leg of a primitive Pythagorean triple is no 
longer than half of the hypothenuse is not particularly restrictive. Moreover, we 
can roughly think of H(n) as a bound on the average orders of a n and b n for each 

n e z >0 . 

On the other hand, we have the following. 
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Lemma 4.3. Let the notation be as above, then C(s) has a pole of order two at 
s = 1 and is analytic for all s G C with SR(s) > 1. Moreover when 5R(s) > 1, C(s) 
has an Euler product type expansion 

^2 



_ I / 4 s -2* + 2 \ -pp _ 1 / 4 s -2* + 2 \ 

1 j Mj ~2l 4- -2- J H — 1 _ 2 V 4 s + 2- J ma 

V 7 prime V / S V , 

where £(s) is </ie Riemann zeta function. 

Proof. Let us consider the Dirichlet series 2£(s), then 

2L(n) 2w(n) Y^ 2 " ( ™ ) Y^ 2 " (n) 1 Y^ 2"( 2 ") 



. OT n\ O u: \ n ) gulnl num) 1 

2£ (s) = = v +y — = y — + 1 y 

n s ^ n s ^ n s ^ n s 4 s ^ n s 

n=l 2\n 4|n 2|n n=l 

2fn \ 2fn 2|n 

/ 2 \ v-^ 2"(") 1 v-^ 2"(") 



4 s / ^ n s 4 s ^ n s 

7 2\n 2\n 

On the other hand, 



2|n n=l \ 2\n 2\n 



and so 



2 " (n) 2 2w(n) 

( > 2v n s - 2 s - 1 1^ ' 



2|n 2fn 
Combining (gSJ) and (gSJ), we obtain 

,50, ^.(f^v; 2 "'"' 

Now define 



4 s - 2 s ^ n s 

2\n 



Li{n) = 



if2|n 
2"(") 2 \ n. 

It is easy to see that Li(l) = 1 and Lj is multiplicative, i.e. if gcd(m,n) = 1 then 
L\{mn) = Li(m)Li{n). Therefore, by Theorem 286 of [TTj . 



- ^-n if ■ n H£<- 

2fn n— 1 / prime \k—0 / 1^2 prime \ fe=l 

< 5i » - n (i^-0= n 

Z^2 prime /^2 prime 

when ffi(s) > 1. Moreover, by Theorem 301 of [11], 

[bA) m^-l^-^r- 11 jrrj- 

n—1 I prime 



k s 
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Now ((47]) follows by combining pj]). ([ST]) , and ((52j) . Moreover, £(s) 2 /C(2s) clearly 
has a pole of order two at s = 1, and is analytic for all seC with 5R(s) > 1. This 
completes the proof. □ 

Remark 4.2. Since £(s) is the sum over all the legs of primitive Pythagorean triples, 
short and long, and it is easy to see that for each (x, y, z) e ty' , max{s, y} > 

Lemma B~B1 combined with Remark |4 . 1 1 imply that Yl( x y z)e<p' m in{a; y} s mus t have 
a pole of order 2 at s = 1 . 



5. Approximating the hexagonal lattice 



In this section we will talk about circle packing density corresponding to simi- 
larity classes of lattices in WR(Z 2 ). Our goal is to prove Theorem ll.6l We do it by 
first proving the following slightly more technical lemma, from which the theorem 
follows easily. 



Lemma 5.1. Let A,B,C be matrices as in Ii26\) . For each k G Z>o, let 

(53) {Pk,t k ,q k ) = (CA) k CB(l, 0, 1) g V, 
then pk > tk, 2|<i,, and C{pk,qk) G C%. Moreover, 

(54) , _ or— 2 Qk-i 

and so 



** = \IiI-pI = 



(55) 



< 



V3 
2 



Pk_ 

qk 



< 



2VS 



0, as k 



oo, 



qk 



(2 + VS)q k (2 + V3)q 2 k 
and more precisely 

(56) 241(7 + 4%/3) fe ~ 1 = 241 x (13. 928. ..) k ~ 1 < q k < 241 x (13.947)* -1 . 
Hence, by I119\) , for each such C{pk,q k ) the corresponding circle packing density is 



(57) 

SO 8m 



1 



12 \ 1 



Pk,Qk 



J Pk,q k ^ 

as in \44\ ) satisfies 



723(7+4 v^)*- 
: 0.9069... = 



4p fe 



< 



1 



12 \ 1 



723x(13.947) fc ~ 1 



qk 



8 (Ah) as k — > oo, and the quadratic form Q Pk , qk (x, y) 

qk - I s 



qk 



xy + y 



(58) — ► Qh(x,y) :— x 2 + xy + y 2 , as k — > oo, 

where Qh(x,y) is the norm form of Ah with respect to the basis matrix as in i20\) . 

Proof. We start by proving (|54"|) . Let {pk,tk, qk) be given by (f53")) . then 

(59) 



4 






1 


3' 




4 
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We argue by induction on k. First notice that p\ — 209, t\ — 120, and q\ = 241, so 
thatpi > ti, 2\ti, and ([54]) is satisfied. Now assume this holds for (pfe-i, ifc-i, ?k-i)- 
By O, 



(60) 
and so 



7pfe_i - 4£ fe _i + 8q fc _i 
4p fe _i -tk-i + 4q fc -i 
3fe_i - 4t fe _i + 9q k -i. 



Pk = t k + (3pfc_i - 3tfe_i + 4g fc _i) > i fc) 
since > t k -i, as well as 

*fc = 4(pfe_i + q fe _i) - tk-l 

is divisible by 2, since 2|tfc_i, and finally 



qk - 1 



4p fe _i - t k -i + Mk-i 



Qk-1 - 1 



t k -: 



tk, 



since — - = tk—i- The conclusion follows by induction 



— 1 _ 

Next we derive (|55|) from ([54]) . Notice that by squaring both sides of ([54]) and 
rearranging terms, we immediately obtain 



(61) 



Pk 

qk 2 _ 



Pk V3 
.» 2 . 



29^ ' 



and since — < 1, we have 

Ik 



Pk 

qk 



V3 

2 



> 



9k 



(2 + V5)^' 

which is the lower bound of (|55|) . For the upper bound, we rewrite (I5T1) as 



ft 



V3 
2 



9fe -2 



< 



1 



< 



qk(2pk + V3qk) 2pk + V3qk 2^3 qk ' 

since V3qk < 2p k - It is also clear that q k — > 00 as fc — > 00. 

To prove ([5"6"|) . we first notice that q% = 241. Moreover, by ([55"]) . the sequence 
Pk/qk is monotone decreasing and converges to \/3/2, therefore 

\/3 < Pfc < Pi = 209 
2 ~ q k ~ qi 241 ' 



for every fc > 1. Then, by ([60]) and ([54]), 

% = 8p fe _i - 4i fc _ a + 9 gfe _i > (7 + 4V% fe _i + 2 > (7 + 4\/% fc _ 1 , 

and 

8 x 209 s 



q k = 8p k - 



4t 



fe-i 



■9« 



fc-i 



< 7 



241 



q fc _i + 2 < 13.947 x q k -x. 



The inequalities l|56[) follow by induction on k. 

To prove (|57[) . notice that upper bound ([53]) implies that 
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where the last inequality is obtained by applying by the lower bound of (|56j) . Then 
the lower bound of (|57| follows. To obtain the upper bound of (|57|) . combine the 
lower bound of (|55p with the upper bound of (|56p in a similar manner. 

Finally notice that (|58|) follows immediately from (|54|) and the fact that qk — > oo 
as fc — > CO, and this completes the proof. □ 

Proof of Theorem ! 1 . 6i Let (A Pfcj9fc ) be the sequence of similarity classes correspond- 
ing to the triples (pk, tk, qk) as defined in (|53|) . then (|55|) guarantees convergence of 
this sequence to the similarity class (Aft) with respect to the metric d s on Sim(]R 2 ), 
and also implies (|2T|) . since qk > qi = 241. The fact that qk = 0(14 fc ) follows 
immediately from (|56|) . To prove (|22p . assume that there exists some similarity 
class (Ap. q ) ^ (Aij) such that 

P 1 
4(A,,A M ) = ---< — 

which implies that (3-\/3p — l) 2 < ^-g 2 , and therefore 

2 2 /1\ 2 27g 2 - 24\/3> + 4 

9 " P ""UJ 27 ' 

where _ 

27g 2 -24^p + 4 2 8p /- 

since if g > 1, then q > 13. Hence 



\Jq 2 -p 2 > 2' 



which contradicts the fact that either (p, -\/g 2 ~ p 2 ,q) or (^/g 2 — p 2 ,p,q) is in T 5 , 
and so (f2"2"|) must be true for each similarity class of the form (A Pi9 ). This completes 
the proof of the theorem, rj 

Finally, Corollary 11.71 follows immediately from (|57|) . 

The approximation result of Theorem 11.61 is also interesting since the similarity 
class (Aft) has a number of important properties: besides providing the optimal 
circle packing and minimizing Epstein zeta function, as mentioned in Section 1, it 
also solves the related minimization problem for the height of flat tori in dimension 
2 (see [S] for details), as well as the quantizer problem in dimension 2 (see jl] for 
details). Let us also recall that a lattice A is called perfect if any real symmetric 
matrix A in the corresponding dimension can be represented as 



\ " t 



A = a x xx 

xeS(A) 

where 5(A) is the set of minimal vectors of A as in Section 1, each x is written 
as a column vector, and each a x is a real number. It is not difficult to see that 
for a lattice A in R 2 to be perfect, the cardinality of S(A) must be six, meaning 
that the only perfect lattices in R 2 come from (Aft). Moreover, (Aft) is strongly 
perfect, meaning that it supports a spherical 5-design: we say that a lattice A in 
R w (and hence its similarity class) supports a spherical t-design for t £ Z>o if for 
every homogeneous polynomial f(x) of degree < t with real coefficients 

( 62 ) / fW dx = wxv\ E /(*)> 
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where 6*"' is the unit sphere in M. N with the canonical measure dx on it, normal- 
ized so that Jqn-i dx = 1. No other similarity class in Sim(R 2 ) supports a spherical 
5-design (or 4-design), and (Ai.i) is the only other similarity class that supports a 
spherical 3-design (or 2-design); such similarity classes are called strongly eutactic 
(clearly, every lattice supports a 1-design) . For detailed information on perfect and 
eutactic lattices see [16], especially chapter 16 for connections to spherical designs. 



6. DlOPHANTINE APPROXIMATION BY QUOTIENTS OF PYTHAGOREAN TRIPLES 



In this section we first prove Theorem 11.81 It follows immediately from the fol- 
lowing direct consequence of a theorem of Hlawka |12j on simultaneous Diophantine 
approximation by quotients of Pythagorean triples, which we state here. 

Then there exist infinitely many 

such that 



Theorem 6.1. Let x 6 (0,1) be a real number. 
Pythagorean triples (p, yj q 2 — p 

(63) x- P - ^ 



< 



Proof of Theorem \1.8[ Recall that 

d s (A, Ap,,) = |sin0(A) - sin6»(A P:9 ) 
and apply Theorem 16. II with x = sin 0(A). 



sin 0(A) - - 
q 



Moreover, we can say that the set {(A Pj9 ) : p/q € S} of similarity classes of WR 
sublattices of 1? is equidistributed in the set Sim(lR 2 ) of similarity classes of all 
WR lattices in R 2 in the following sense. It is a well known fact that the map 



t 



1-t 2 



t 



is a bijection from the set of rational numbers onto the set of all rational points on 
the unit circle. Ordering Q as the set of Farey fractions induces an ordering on the 
set of rational points on the unit circle, and hence on the set S of y-coordinates of 



such points that fall in the interval 



V^5 



Now, it is a well known fact that Farey 



fractions are uniformly distributed (modi). 

As a side remark, we can also use Theorem 16. II to approximate points on a unit 
circle with rational points on the same circle. 

Corollary 6.2. Let (x, y) be a point on the unit circle. Then either x, y S {0, ±1}, 
or there exist infinitely many rational points (p/q,r/q) on the same circle such that 

(64) max I x — — , y !> X 



Q 



Proof. First notice that it suffices to prove the statement of this corollary for the 
case < x, y < 1, namely the case when the point in question lies in the first 
quadrant, since any other point on the circle can be obtained from those in the first 
quadrant by a rational rotation. Let c be an arbitrary real number in the interval 
(0, 1), then either 



(65) 



< x < \J\ - c 2 < 1, c < y < 1, 



26 



LENNY FUKSHANSKY 



(66) < y < v 7 ! -c 2 < 1, c<x< 1. 

First assume that (|65|) holds. By Theorem l6.11 there exist infinitely many Pythagorean 
triples (p, r, q) with r = \/ q 2 — p 2 which satisfy (|63p . Then: 

2y^ 



> 



(67) 







a; 














1 — — 5 



1 



r 

y- - 

9 



> 



c 1 + 



n+l 



2«/l 



r 

y - - 

9 



The last inequality is true because ^= a ^ s an increasing function in both 

variables for < z, w < 1; since y > c, we can pick g large enough so that r/g 
would have to be sufficiently close to y so that r/q> for some n € Z>o, then 

and v/T^ + < 2 A /l - (^c 2 . Then flUD 



r/g + y > c (l + ^ 
implies: 

(68) 



< 



1 - 



4-y/2 



c 1 



n+l 



9 



Since our choice of c G (0,1) and positive integer n was arbitrary, we can for instance 
choose 



(69) 



2n + 2 



4n + 1 



and take n = 2, in which case, combining (|63|) . (|68|) . and (|69|) . we obtain (|64|) . 

If, on the other hand, holds instead of (JHSJ), simply repeat the above argu- 
ment interchanging x with y and p/q with r/g. This completes the proof. □ 

A related result has also been obtained by Kopetzky in [T5] (also see |14|). 
however his bounds are different in flavor in the sense that the constants in the 
upper bounds depend on x and y. Notice that the bound of Corollary 16.21 can be 
easily extended to any rational ellipse. 

Corollary 6.3. Let (x,y) be a point on the ellipse E, given by the equation 



(D 2 +(! 



where a,b are positive rational numbers. Then either (x,y) = (±a, 0), (0, ±6) 7 or 
there exist infinitely many rational points {p/q,r/q) on the same ellipse such that 



(70) 



< 



2\/2 max{a, b} 



q 



Proof. Notice that the map (x,y) <— > (x/a,y/b) is a bijection between E and the 
unit circle, which takes rational points to rational points. Now apply Corollary 16. 21 
to points of the form (x/a, y/b). □ 
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